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Abstract 


The warming of the air and the earth before the occurrence of earthquakes is a precursor that has 
been known since ancient times. Today, the truth of this myth, with the scientific title of thermal 
anomaly, has been proven in many earthquakes. The present study aims to investigate the occurrence 
of the anomaly in the long-term period of three years before the 2017 Azgeleh earthquake, using the 
simple method of deviation from the moving average. The results showed that one or two weeks 
before the occurrence of the earthquake, a rare or unprecedented daily anomaly occurred in both the 
meteorological parameters of air temperature, soil temperature, and air humidity at the synoptic 
station of Tazeh-Abad, as well as in the satellite parameter of the land surface temperature. These 
thermal anomalies can be related to the occurrence of the earthquake. The spatial investigation of 
the land surface temperature anomaly that occurred 15 days before the earthquake, revealed that on 
that day, the land boundary between the Arabian and Eurasian plates in the west of Iran was warm, 
and a hot area appeared in the center of Iraq with elongation toward the epicenter of the earthquake. 
Other tests showed that this simple anomaly detection method is also able to identify the thermal 
anomaly before the 2003 Bam and 2010 Mohammad-Abad Rigan earthquakes, both in Kerman 
province. The Azgeleh earthquake was a typical earthquake of the active and densely populated 
Zagros zone. This type of research hopes to better understand the precursors of earthquakes in each 
zone, with the aim of realizing the dream of earthquake prediction. 
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حکیده 


گرم شدن هوا و زمین پیش از رخداد زلزله‌هاء پیش‌نشانگری است که از دوران باستان شناخته شده 
است. امروزه درستی این افسانه. با عنوان علمی ناهنجاری حرارتی. در زلزله‌های بسیاری به اثبات 
رسیده است. هدف این پژوهش, بررسی بروز این ناهنجاری در 99 )8 بلندمدت سه MUS‏ پیش از 335 
AS 5I Wat‏ کرمانشاه با روش dala‏ انحراف از میانگین متحرک است. نتایج نشان داد که در یکی دو 
هفته پیش از رخداد این زلزله» هم در پارامترهای هواشناسی دمای هو دمای خاک و رطوبت هوا در 
ایستگاه همدید تازه‌آباد و هم در پارامتر ماهواره‌ای Glee‏ سطح زمین, ناهنجاری روزانۀ کم‌سابقه یا 
بی‌سابقه‌ای رخ داده است. این ناهنجاری‌ها حرارتی می‌تواند با رخداد زلزله در پیوند باشد. بررسی 
مکانی ناهنجاری در دمای سطح زمین که ۱۵ روز پیش از زلزله رخ داده بود. روشن ساخت که در آن 
روز» مرز خشکی دو صفحة عربستان و اوراسیا در غرب ایران. گرم بوده و ناحیه‌ای plo‏ در مرکز عراق 
با کشیدگی به سمت رومرکز زلزله پدیدار شده است. دیگر آزمون‌ها نشان داد این روش dale‏ کشف 
ناهنجاری» قادر به شناسایی ناهنجاری حرارتی پیش از زلزله‌های ۱۳۸۲ بم و ۱۳۸۹ محمدآباد ریگان, 
هر دو در استان کرمان نیز هست. زلزلة ازگله» شاخحصی برای زلزله‌های iiku‏ فعال و پرجمعیت 
زاگرس است. امید این دست پژوهش‌ها. شناخت بهتر پیش‌نشانگرهای زلزله در هر منطقه با هدف 
تحقق رژیای پیش‌بینی زلزله است. 

کلیدواژه‌ها: ناهنجاری حرارتی» پیش‌بینی زلزله» پیش‌نشانگرهای زلزله» clo‏ سطح زمین» روش 


انحراف از میانگین متحرک» زاگرس ایران 
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مقدمه 


زلزله از ویرانگرترین مخاطرات است. بشر دریافته که پیش از زلزله‌های بزرگ» نشانه‌هایی پدیدار می‌شود 
(جدول ۱). هدف از بررسی این پدیده‌ها؛ پیش‌بینی زلزله‌های tel‏ است. امرون هرگونه تغییر 
قابل‌اندازه‌گیری در میانگین مشاهدات محیطی. پیش از رخداد زلزله را با نام پیش‌نشانگر زلزله می‌شناسیم. 
یکی از آن‌هاء ناهنجاری حرارتی است که مهم‌ترین پارامتر آن» دمای سطح زمین LST)‏ است و عموماً یکی 


42 هفته پیش از رخداد زلزله‌های بزرگ» a> > A>‏ سلسیوس در رومرکز افزایش می‌یابد (Zhao, Pan, Sun,‏ 
.Guo, Zhang, & Feng, 2021)‏ 


با آغاز عصر فضا و پرتاب ماهواره‌های دیدبانی زمین در des‏ دوم Sule‏ بیستم. علم سنجش‌ازدور پا گرفت 
که هنر و فن‌آوری کسب اطلاعات از پدیده‌ها از راه دور است. در ده ۱۹۸۰ برای نخستین ob‏ تغییر دمای 
سطح زمین پیش از زلزله‌های آسیای میانه در تصاویر حرارتی دیده شد؛ و این‌چنین تاریخچۀ کاربرد 
سنجش ازدور حرارتی در بررسی پیش‌نشانگر ناهنجاری حرارتی در زلزله‌ها گشوده شد(2010 (Tronin,‏ 

کلید پیش‌بینی زلزله» پایش پیش‌نشانگرها و کشف ناهنجاری در آن‌ها با تحلیل سری زمانی است. هر Am‏ 
شدت اهنجاری بیشتر wah‏ زلزله‌ای با بزرگای بیشتر و يا عمقی کمتر رخ می‌دهد (Shen, Zhang, Hong, Jing,‏ 
Zhao, 2013)‏ & اما منتقدان معتقدند مبنای نظری پیش‌نشانگرها و مکانیسم ناهنجاری‌ها کاملاً روشن نیست؛ و 
زلزله‌های بزرگ نادرند و قابل تکرار نیستند. همچنین داده‌های ماهواره‌ای و الگوریتم‌های سنجش‌ازدوری دارای 
عدم قطعیت‌اند و با Gla le‏ فراوانی روبرو هستند(2021 .(Picozza, Conti, & Sotgiu,‏ 


حدول ۱- فهرست AS‏ پیش‌نشانگرهای رایج زلزله(2021 (Conti, Picozza, & Sotgiu,‏ 


Tablel- General list of common earthquake precursors 


P‏ - وقوع پیش‌لرزه و تغییر لرزه‌حیزی - تغییر دما و رطوبت سطح و اعماق زمین 
Se:‏ + 
- جابجایی و دگرشکلی زمین - تغییر رسانایی و مقاومت خاک 
p‏ - تغییر دبی رود و چشمه - pad‏ دمای col‏ چاه رود. چشمه و bos‏ 
آب کره ۱ : 
- تغییر سطح اب‌های زیرزمینی - تغییر عناصر شیمیایی اب چشمه و چاه 
- تغییر دما و رطوبت و فشار هوا - تغییر محتوای هواویز 
هواکره 
- تغییر محتوای گازهای جو - تشکیل ابر زلزله 
1 - تغییر محتوای کل الکترون و شب‌تابی - آشفتگی میدان الکتریکی پا مغناطیسی زمین 
A‏ ۲ 
= تغییر گرانش زمین - تغییر رفتار جانوران 


1 Land Surface Temperature 
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گرم شدن هوا و زمین پیش از زلزله‌هاه از دوران چين باستان شناخته شده بود )2010 (Feng & wu‏ 
اعتقاد بر این است که تنش موجود در پوستۀ زمین و افزایش OF‏ پیش از زلزله. سبب افزایش LST‏ و درنتیجه 
افزایش تابش مادون‌قرمز حرارتی TIR)‏ می‌شو د )2009 (Saraf, Rawat, Choudhury, Dasgupta, Das,‏ 
این افزايش تابش TIR.‏ که bug‏ ماهواره‌ها قابل تشخیص بوده و بروز Ol‏ پیش از زلزله‌ها بارها گزارش شده 
است. ناهنجاری حرارتی خوانده می‌شود(2005 (Freund, et al.,‏ 

تاکنون چندین فرضیه برای توجیه گرمای پیش از زلزله ارائه شده است؛ ازجمله تراوش گاز و بخارآب از 
زمین در اثر باز شدن حفره و شکاف‌های درون زمین و انتقال گرمای همرفتی با تشدید اثر گلخانه‌ای محلی 
(Qiang, Xu, & Dian, 1991)‏ ایجاد LF‏ در اثر اصطکاک بین بلوک‌های 145 )2017 (Bhardwaj, et al.,‏ « 
گرمازایی در واکنش یونیزاسیون هوا توسط گاز مادون منتشرشده از زمین )2004 (Pulinets,‏ و جدا شدن 
الکترون‌های آنیون اکسیژن از سنگ‌های آذرین که سبب تابش امواج مادون‌قرمز میانی MIR)‏ و ثبت آن در 
تصاویر حرارتی به‌صورت LŚ‏ می‌شود(2002 (Freund,‏ 

بررسی ناهنجاری حرارتی در زلزله‌هاء ادبیاتی غنی دارد )2015 (Tramutoli, et al.,‏ نخستین مورد در 
تاریخ ple‏ متعلق به )1988 (Gornyi, Salman, Tronin, & Shilin,‏ است که بروز آن پیش از زلزله‌های 
آسیای میانه را در تصاویر حرارتی ماهوارۀ NOAA‏ نشان دادند.(1997 (Qiang, Xu, & Dian,‏ با تصاویر 
ماهوارة Meteosat‏ به افزایش دمای زمین VU‏ درجۀ سلسیوس در هفتۀ منتهی به زلزله‌ای در ce oem‏ بردند. ) 
Tronin, Hayakawa, & Molchanov, 2002)‏ با بررسی هفت سال تصاویر Sole‏ ماهوارءٌ NOAA‏ در 
شرق del‏ بین ناهنجاری‌های مثبت ule‏ در گسل‌ها و رخداد زلزله‌های کم‌عمق در چین و ژاپن؛ 
همبستگی خوبی يافتند. 

(Ouzounov, & Freund, 2004)‏ برای اولین بار با محصولات MODIS ix‏ ثابت کردند که V‏ روز 
پیش از 335 گجرات هند. LST‏ در رومرکز تا ۶ de yo‏ سلسیوس افزایش aab‏ است. )2008 (Saraf, et al.,‏ 
با کاوش ده زلزلۀ ایران ازجمله زلزلٌ بم روی تصاویر حرارتی ماهوارة NOAA‏ دریافتند که ۱۰ تا ۱ روز 
پیش‌ازاین زلزله‌ها؛ افزایشی در LST‏ رومرکز با شدت Y‏ تا ۱۳ درجۀ سلسیوس رخ داده است. 

(Saradjian, & Akhoondzadeh, 2011)‏ با محصول LST‏ سنجندة MODIS‏ بروز ناهنجاری حرارتی. 
۰ تا ۱ روز پیش از زلزله‌های بم» زرند و بروجرد را با شدتی تا Y‏ درجه سلسیوس شناسایی کردند. 
(Akhoondzadeh, 2014)‏ بروز ناهنجاری LST‏ در ^ روز پیش از 4335 سراوان را کشف کرد. (Bellaoui,‏ 
Hassini, & Bouchouicha, 2017)‏ با بررسی پنج سال تصاویر حرارتی شبگاه سنجندة MODIS‏ و آمار 


! Thermal Infrared 
? Mid-Infrared 
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دمای هوای سه ایستگاه هواشناسی نزدیک» ناهنجاری حرارتی مثبتی را به‌طور هم‌زمان در هر دو سری داده 
پیش از زلزله‌ای در الجزایر یافتند. 

(Barkat, et al., 2018)‏ در محصول LST‏ سنجندة MODIS‏ افزایشی ۲ تا ۱۰ درجۂ سلسیوس در یکی 
دو هفته پیش از سه زلزلة پاکستان یافتند.(2019 Khalili, Alavi Panah, Abdollahi Eskandar,‏ ) با کاوش 
همین محصول در 8599 ده‌ساله. رومرکز زلزله سراوان را در سال رخداد این زلزله و در روزهای پیش و 
پس‌ازآن» ناهنجار دیدند. 

(Genzano, Filizzola, Hattori, Pergola, & Tramutoli, 2021)‏ با بررسی ۱ سال تصاویر حرارتی 
شبگاه ماهواره‌های MTSAT‏ دریافتند که 1۲ درصد ناهنجاری‌های حرارتی رخ‌داده در ژاپن» به زلزله‌ای با 
بزرگای دست‌کم 1 ختم شده است. عین همین پژوهش را )2022 (Filizzola, et al,‏ با بررسی ۱۲ سال 
تصاویر حرارتی شبگاه ماهوارۀ Meteosat‏ روی ترکیه انجام دادند و به نسبت VE‏ درصد برای زلزله‌های با 
بزرگای دست‌کم ۶ رسیدند. 

(Askari, Hafezi, Rahimi tabar, & Ansari, 2010)‏ در تصاویر حرارتی ماهوار؟ NOAA‏ افزایش 
شدیدی در 1 روز پیش و در رومرکز 335 راور دیدند. (Choubsaz, Akhoondzadeh, & Saradjian,‏ 
)2015 نیز با همین تصاویر بروز ناهنجاری حرارتی در یکی دو هفته پیش از سه زلزلة اهر SIS‏ و سراوان 
را یافتند. )2017 (Saber Mahani, & Sepahvand,‏ با محصول LST‏ شبگاه MODIS ài‏ بروز 
ناهنجاری حرارتی در Y‏ روز پیش از 335 محمدآباد ریگان را شناسایی کردند. (Khoshgoftar, Saradjian,‏ 
)2021 در محصول LST‏ سنجندةٌ MODIS‏ و سه پیش‌نشانگر ORs‏ انواع ناهنجاری در دو U5‏ بم و ازگله 
را تشخیص دادند؛ چنان که $36 زمانی و بزرگای هر دو قابل پیش‌بینی بوده است. 

oly!‏ در کمربند 335 33 آلپ هیمالیا و در برخوردگاه صفحة عربستان با اوراسیا قرار دارد و زلزله‌خیز است. 
واپسین زلزلة شدید ایران» زلزلة ۱۳۹۹ A 3I‏ کرمانشاه است. در این پژوهش. یک بار دیگر بروز پیش‌نشانگر 
حرارتی در این زلزله. ولی برای اولین بار در دوره‌ای بلندمدت و با روشی ساده و با کمک داده‌های ایستگاه 
همدید. بررسی خواهد شد. زلزلة ازگله. شاحصی برای زلزله‌های منطقة لرزه‌ای زاگرس است. هر دو سال 
یک بار» زلزله‌ای با بزرگای دست کم V‏ سبب جان باختن و بی‌خانمانی ساکنان این منطقه می‌شود. منطقة 
زاگرس» ده استان و یک‌چهارم جمعیت کشور را در خود جای داده است. هر شناختی از فرآیند 
پیش‌نشانگری زلزله‌ها با رویکرد هشدار یا پیش‌بینی, می‌تواند منجر به حفظ جان و مال مردم و مدیریت بهتر 


بحران گردد. 


سال سیزدهم حیدری و همکاران» بررسی pirum‏ ناهنجاری حرارتی در زمین‌لرزه ۱۳۹۹ از aS‏ کرمانشاه av‏ 
مواد و روش‌ها 
منطقة موردمطالعه 


مرکز ثقل رومرکزهای زلزله در گزارش‌ها (جدول dY‏ به‌عنوان مرکز منطقۀ موردمطالعه در نظر گرفته شد. 
این نقطه در ۲ کیلومتری جهت ENE‏ شهر ازگله قرار دارد. مرز منطقه به اندازۀ یک انحراف معیار از 


پراکند گی )29 La 5S‏ دورتر از رومرکزهای پیرامونی oA AS‏ شد. منطقة موردمطالعه. ۷ کیلومتر طول و 
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شکل Y‏ — موقعیت منطقة موردمطالعه. رومرکزهای 5 از گله و ایستگاه سینوپتیک تازه‌آباد. 


Fig.1. The location of the study area, epicenters of Azgeleh earthquake 0 0 
synoptic station 


زلزلة موردمطالعه 


روز یکشنبه» ۱۲ نوامبر ۲۰۱۷ ساعت ۱۸:۱۸ به‌وقت جهانیء برابر با ۲۱ آبان ۱۳۹۲ ساعت ۲۱:۶۸ به‌وقت 


کرمانشاه و در چند کیلومتری مرز عراق به وقوع پیوست. 


بیش از ٩‏ هزار نفر مجروح شدند. همچنین. جند شهر به‌طور کامل یا 


ol‏ دست کم Wwe‏ نفر Oe‏ باختند و د 


۹۸ نشریه جغرافیا و مخاطرات محبطی» جلد OY‏ شماره ۲ تابستان ۱۶۰۳ 


جزئی تخریب شد و ۰ هزار واحد مسکونی روستایی TAM]‏ حدی دید و در مجموع 15.55 ۶۳۰ هزار نفر 
تحت yb‏ قرار گرفت (2021 (Maghsoudi & Moshtari,‏ و در پی عدم تاب‌آوری مردم. پیامدهای ناگوار 
اجتماعی. خانوادگی. اقتصادی. بهداشتی و روانی به بار آورد(2021 (Shahbazi, Papzan, & Gholami,‏ 


ویژگی‌های لرزه‌شناسی این زلزله از آخرین نسخهٌ گزارش‌های بین‌المللی گرفته شد (جدول ۲ و شکل A‏ 


جدول ۲- ویژگی‌های لرزه‌شناسی زلزلة ۱۳۹ ازگله. 


Table2- Seismological characteristics of the 2017 Azgeleh earthquake 


زمان طول عرض عمق بزرگا 
(M) (km) N) (E) (UTC)‏ 


WK ۱ ۳٤/۷۷ 0/۷17 ۸۷4 IRSC‏ مرکز لرزه‌نگاری oy gS‏ دانشگاه تهران 


منبع گزارش 


TIAN EOE ۱۸۱۸۱۸۸ HEES‏ ۰ ۱۸ ۳ پژوهشگاه بین‌المللی زلزله‌شناسی ایران 
USGS‏ ۱۸۷/۲ 10/۹۰ ۳/۹۱ ۰ ۱۹ ۷/۳ سازمان زمین‌شناسی آمریکا 

Yé ۳ ۳ EMSC‏ ۳ مرکز زلزله‌شناسی ارویا و مدیترانه 
۳٤۸۵ ۵/۹۲ ۷/۲ GFZ‏ ۲۲ ۷/۳ مرکز پژوهشی علوم زمین آلمان 


۱ VY ۲۰ ۳/۸۶ ۵/۸۷ ۱۸:۱۸:۷۲  نیگنایم‎ 


دریافت داده‌های هواشناسی 


نزدیک‌ترین ایستگاه هواشناسی همدید. در ۱۱ کیلومتری منطقة موردمطالعه» ایستگاه تازه‌آباد با کد 
Lots KEK‏ ۹۹۶۲۸ و ارتفاع ۱۲۲۹ متر از سطح Uo»‏ است (شکل ۱). انواع پارامترهای sles‏ هوا 
e (AT)‏ رطوبت هوا و دمای خاک این ایستگاه برای دورۀ سه Ble‏ ۲۰۱۵ تا پایان ۲۰۱۷ از سازمان هواشناسی 
کشور IRIMO)‏ دریافت شد. داده‌های این ایستگاه مفقودی کمی دارد ولی از لحظة زلزله تا دو هفته قطع 


می‌شود. 


دریافت داده‌های ماهواره‌ای 


دو "MODIS iu‏ در سال‌های ۱۹۹۹ و ۲۰۰۲ Lug‏ سازمان ملی هوانوردی و فضایی آمریکا 
NASA)‏ روی ماهواره‌های Terra‏ و Aqua‏ با هدف GEL‏ تغییرات زمین در مدار قرار گرفتند. آن‌ها در 
ارتفاع ۷۰۵ کیلومتری و در مداری شبه‌قطبی و خورشیدآهنگ با زمان بازدید VV‏ روزه قرار دارند و به ترتیب 
Atmospheric Temperature‏ ! 

? [slamic Republic of Iran Meteorological Organization 


? Moderate Resolution Imaging Spectroradiometer 
^ National Aeronautics and Space Administration 


سال سیزدهم حیدری 3 همکاران بررسی تش نشانکر ناهنجاری حرارتی در زمین‌لرزه ۱۳۹۹ از aS‏ کرمانشاه ۹۹ 


در مدار نزولی با گذر ۱۰:۳۰ و مدار صعودی با گذر ۱۳۳۰ به‌وقت محلی بر استواء با توان تفکیک مکانی 
دست کم ۱ کیلومتر» در YA‏ باند تصویر برمی‌دارند(2002 (Justice, et al.,‏ 

محصول دما و گسیلمندی سطح زمین در ماهوارة Aqua‏ با نام MYDIIAI‏ به‌طور روزانه با توان 
تفکیک مکانی ۱ کبلومتر و در گسترة جهانی تولید می‌شود. در این پژوهش از LY‏ دمای سطح زمین (LST)‏ 
شبگاه از نسخهٌ 7 استفاده شد. اگرچه این محصول توان اندازه‌گیری در آسمان ابری را ندارد؛ ولی W‏ درصد 
از شب‌هاء ابرناکی منطقۀ موردمطالعه کمتر از ۱۰ درصد است (شکل ۲). این محصول در انواع پوشش‌های 
زمین و شرایط مختلف آب‌وهوایی با داده‌های زمینی اعتبارسنجی شده و صحت LST‏ شبگاه ol‏ با خطای 


جذر میانگین مربعات (RMSE)‏ کمتر از Y‏ کلوین به تأیید رسیده است )2019 (Duan, et al.,‏ 


A: Daily Nighttime Cloudiness Percentage 
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شکل ۲- وضعیت ابرناکی شبگاه منطقۀ موردمطالعه. -A‏ روزانه؛ -B‏ هیستوگرام اپرناکی؛ -C‏ میانگین ماهانه. 


Fig.2. The nighttime cloudiness condition of the study area: A- Daily; B- Cloudiness 
histogram; C- Monthly average 


روش انجام پژوهش 

از همان نخستین پژوهش‌ها در این زمینه. پژوهشگران به برتری تصاویر حرارتی شبگاه نسبت به روزگاه 
آگاه بوده‌اند. تصاویر حرارتی شبگاه وابستگی کمتری به گرمای تابشی خورشید دارند و GES‏ درونی زمین 
را بهتر OUS‏ می‌دهند. همچنین در روز. بخشی از امواج حرارتی خورشید در اثر پخش یا انعکاس به سنجنده 
می‌رسد که خطاساز است(2005 (Saraf & Choudhury,‏ آن‌ها که هر دو سری را در نظر گرفته‌اند» دیده‌اند 
که ناهنجاری حرارتی نخست در سری شبگاه „alb‏ شده(2008 (Saraf, et al.,‏ و Ol GAS‏ در سری روزگاه 
خفیف‌تر است؛ چنان که گاه تنها در سری شبگاه مشاهده می‌شود & (Mohamed, Gahalaut, Sekertekin,‏ 


Inyurt, 2021) 


! Root Mean Square Error 


نشریه Lal ae‏ و مخاطرات محیطی. جلد OY!‏ شماره ۲. تابستان ۱۶۰۳ 


بنابراین از دیرهنگام‌ترین تصاویر LST‏ که برای باند شبگاه محصول 1۷۷۲11۸1 از ماهوارة Aqua‏ و 
برای حوالی زمان محلی Y‏ بامداد است؛ استفاده گردید. به دلیل استفاده از محصول LST i36‏ مراحل 
پیش پردازش تنها در ۶ مرحلهٌ ساده در سامانۂ موتور زمین گوگل ples! (GEE)‏ شد (شکل ۳): 

۱- فراخوانی محصول: سری محصول )9 LST lj‏ در ijb‏ زمانی مورد نظر فراخوانی شد. 

۲- برش تصاویر: همه تصاویر طبق مرز منطقة موردمطالعه برش خورد. 

۳- حذف تصاویر ابری: تصاویر با آستانۂ ابرناکی ۵۰ درصد کنار گذاشته شد. 


-t‏ میانگین گیری: برای استخراج LST‏ روزانه از هر تصوير میانگین گرفته شد. 


1- Call Daily LST Series 2- Clip to the Study Area 3- Remove Too Cloudy LSTs 4- Average Every LST 


Fig.3. Preprocessing steps of LST product in order to detect thermal anomaly 


برای اعتبارسنجی محصول LST‏ شبگاه از داده‌های ایستگاه هواشناسی همدید تازه‌آباد استفاده شد. بازه 
زمانی سه JUS‏ موردمطالعه یعنی از آغاز ۲۰۱۵ تا پایان ۲۰۱۷ شامل ۱۰۹۷ روز می‌شود. در این «mde‏ 
ایستگاه هواشناسی ۱۲۰ روز مفقودی و محصول LST‏ نیز ۳۶۷ روز SU pl‏ در OF‏ پیکسل داشت. این دو 
۰ روز برابر با M‏ درصد 556 زمانی سه‌ساله را دادۀ هم‌زمان دارند. 

ازآنجایی که محصول LST ilja‏ از تصاویر برداشتی در گذرهای مختلف تولید می‌شود. محصولات. 
زمان محلی متفاوتی بین ۱:۳۷ تا ۳:۱۶ با متوسط ۲:۲۳ و انحراف معیار aids YA‏ داشت. چون ایستگاه 
تازه‌آباد. دمای سطح زمین را تنها در صبح‌ها برداشت کرده است. برای اعتبارسنجی محصول LST‏ از دمای 
هوای میان‌یابی‌شده در همان زمان تصویربرداری استفاده گردید و همبستگی آن‌ها با ضریب تعیین ۰/٩۱‏ تأیید 


شد. در اکثر شب‌هاء cleo‏ هوا به‌طور متوسط ۶/۳ درجۀ سلسیوس گرم‌تر از دمای سطح زمین است (شکل 
Ak‏ 


! Google Earth Engine 


سال سیزدهم حیدری و همکاران» بررسی پیش‌نشانگر ناهنجاری حرارتی در 5 Jose‏ $5 ۱۳۹۲ ازگلۀ کرمانشاه ۱۰۱ 
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شکل ۶- اعتبارسنجی محصول LST‏ در سنجنده MODIS‏ با دمای هوا در ایستگاه تازه آباد. 


دب 
a‏ 


Fig.4. Validation of MODIS LST product with air temperature at Tazehabad station 


در شب. زمین حرارت خود را تابش کرده و خنک‌تر از هوا می‌شود. سایر پژوهش‌ها نیز همبستگی قوی 
و مثبت این دو پارامتر را با رابطه‌ای حطی نشان داده‌اند(1992 (Bevis, et al.,‏ که البته شیب ol‏ متأثر از طبقة 


پوشش اراضی و به‌ویژه نوع پوشش گیاهی منطقه است(2011 (Mildrexler, Zhao, & Running,‏ 
نتایج و c‏ 


داده‌های هواشناسی ایستگاه همد 

سری داده‌های پارامترهای هواشناسی ایستگاه تازه‌آباد به‌عنوان شاهدی از منطقه. در سه بخش بررسی شد. 
در مورد هر پارامت ناهنجاری روزانه به‌صورت میزان انحراف از میانگین متوسط M‏ روزه محاسبه گردید. 
دمای هوا 


das cle‏ روزانه در ۸ پارامتر سه‌ساعته و ۲ پارامتر clea‏ کمینه و دمای بيشینه تبت می‌شود. در بیشتر این 
پارامترها و به‌ویژه در پارامترهای شبگاه» ناهنجاری مثبتی پیش از زلزله دیده می‌شود که چشمگیرترین ol‏ در 
elk‏ دمای Ana S‏ رخ داده است (شکل 0( 


۱۰۲ نشریه Lal ae‏ و مخاطرات محیطی. جلد VY‏ شماره ۲. تابستان ۱۶۰۳ 


A: TazehAbad Minimum Air Temperature 
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شکل 0— سری داده دمای کمینه در ایستگاه تازه‌آباد: -A‏ داده‌های روزانه؛ -B‏ ناهنجاری روزانه. 


Fig.5. Minimum temperature data series at Tazehabad station: A- Daily data; B- Daily 
Anomaly 


slays‏ بالا ھان dis ce‏ که در این دوا all utes‏ کین atas aos ed ad‏ درج مامیرس 
برای YA‏ اکتبر ۲۰۱۷ یعنی ۱۶ روز پیش از زلزله بوده است. )2022 (Jing & Singh,‏ نیز در داده‌های سنجندة 
5 دای هوا در ۲۷ اکتبر را در رومرکز و tte‏ بخش‌های خاورمیانه دارای ناهنجاری مثبت یافته‌اند. 

همچنین بررسی داد ایستگاه‌های هواشناسی در زلزلۀ قوی مکزیک نشان داده که یکی دو هفته پیش 
از زلزله‌هاء ناهنجاری مثبت در دمای هوا و دامنۀ دمای روزانه € داده است (Dunajecka & Pulinets,‏ 
)2005 بررسی مدل هواشناسی در زلزله‌ای در چین نیز ناهنجاری مثبت دمای هوا در ۱۲ روز پیش از زلزله را 
نشان می‌دهد که این گرما از زمین سرچشمه گرفته و سپس به لایه‌های SVL‏ جو ume‏ شده و سرانجام 
ناپدید شده است(2021 (Zhang, Meng, Wang, Lu, Hu,‏ بررسی چهارده 40305 قوی آسیای میانه نیز OUS‏ 
داده که در dew‏ موارد. cles‏ هوا در ماه رخداد زلزله از dos‏ ماه‌های مشابه در دوره‌ای چند ده‌ساله گرم‌تر بوده 
است )1986 (Milkis,‏ 


دمای خاک 


دمای خاک روزانه در ۱ پارامتر برای سطح زمین در ساعت V‏ جهانی و در ۱ پارامتر برای اعماق ۵ cha‏ 
۰ ۳۰ ۵۰ و ۱۰۰ سانتی‌متر در ساعت ٩‏ جهانی ثبت می‌شود. با اینکه هیچ کدام از این پارامترها برای شبگاه 
ALL LEM LN:‏ اک( A‏ ناهنجاری 


مثبت کم‌سابقه‌ای پیش پیش از زلزله دیده می‌شود. 


! Atmospheric Infrared Sounder 


A: TazehAbad Soil Surface Temperature at 06:00 UTC 
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شکل 7 سری cles dolo‏ سطح خاک روزگاه در ایستگاه تازه‌آباد: -A‏ داده‌های روزانه؛ -B‏ ناهنجاری روزانه. 


Fig.6. Daytime soil surface temperature data series at Tazehabad station: A- Daily data; 
B- Daily anomaly. 


طبق نمودار VL‏ در $555 سه AU.‏ پیش از زلزله Ota‏ انحراف بزرگ دمای سطح خاک با ۷/۷ درجة 
سلسیوس برای ۰ odas ۷ ysl‏ ۱۳ روز پیش از زلزله بوده Aud‏ بررسی این پارامتر در ایستگاهی در 


رومرکز زلزله‌ای در مکزیک نیز OLAS‏ داده که میانگین دمای سطح خاک ماهانه در ol‏ ماه از ماه‌های مشابه در 


۰ سال اخیرش گرم‌تر بوده است(2006 (Pulinets, et al.,‏ 


رطوبت هوا 


رطوبت dg‏ روزانه در پارامترهای سه ساعتۀ رطوبت نسبی هوا و فشار بخار ثبت می‌شود. در بیشتر این 
پارامترها و به‌ویژه در فشار بخار )59 ol‏ ناهنجاری جشمگیر مثبتی پیش از زلزله دیده می‌شود؛ که اوج آن در 
پارامتر فشار بخار در ساعت ٩‏ جهانی (شکل (v‏ نمود دارد. 


۱۶۰۳ شماره ۲. تابستان‎ OY! و مخاطرات محیطی. جلد‎ Lal ae نشریه‎ yt 


A: TazehAbad Vapor Pressure at 09:00 UTC 
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شکل -N‏ سری ols‏ فشار بخار ساعت ۹ حهانی در ایستگاه تازه آباد: -A‏ داده‌های روزانه؛ -B‏ ناهنحاری روزانه. 


Deviation (hpa) 


Fig.7. Vapor pressure data series at 9 UTC at Tazehabad station: A- Daily data; B- Daily 
anomaly. 


GL‏ نمودار AL‏ بیشترین انحراف فشار بخار با ۷/۱ هکتوپاسکال برای ۵ نوامبر ۲۰۱۷ یعنی ۷ روز پیش از 
زلزله بوده است. )2019 (Akhoondzadeh, De Santis, Marchetti, Piscini, Jin,‏ نیز در مدل جهانی جو 
۲ محتوای بخار آب را ٩‏ روز پیش از این زلزله ناهنجار یافته‌اند. افزايش ناگهانی این پارامتره ۱۱ 
روز پیش از زلزله‌ای در چین و پس از jay‏ ناهنجاری حرارتی در آن گسل )2009 (Liu, Cui, Wu, Wang,‏ 
7 روز پیش از زلزله‌ای در آلاسکا و ۱ تا ۱۶ روز پیش از ۶۰ Jj‏ ساحلی نیز رخ داده (Singh, Gul‏ 
Cervone, Singh, & Kafatos, 2007)‏ 


کشف ناهنجاری در دمای سطح زمین شبگاه 


ناهنجاری زمانی 


ناهنجاری زمانی در داده‌های روزانه» پس از کنار گذاشتن تصاویر ابری با آستانة ابرناکی ۵۰ درصد. 
به‌صورت میزان انحراف از میانگین متحرک ۲۱ روزه محاسبه گردید (شکل (A‏ حذف نکردن تصاویر ابری از 
سری زمانی می‌تواند موجب سوگیری در LST‏ و نمود ناهنجاری GIS‏ شود & (Blackett, Wooster,‏ 
Malamud, 2011)‏ 


! European Centre for Medium-range Weather Forecasts 


سال سیزدهم حیدری و همکاران» بررسی پیش‌نشانگر ناهنجاری حرارتی در زمین‌لرزۀ ۱۳۹۲ ازگلۀ کرمانشاه ۱.0 


A: Daily Nighttime LST Averages in Region A (Cloudiness <50%) 
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شکل ۸- سری زمانی LST‏ شبگاه در منطقة موردمطالعه: -A‏ داده‌های روزانه؛ -B‏ ناهنحاری روزانه. 


Fig.8. Nighttime LST time series in the study area: A- Daily data; B- Daily anomaly 


Gb‏ نمودار VL‏ در این دور سه JUS‏ پیش از زلزله» بیشترین انحراف LST‏ شبگاه در منطقة موردمطالعه 
با WY‏ درجۀ سلسیوس و نمرة استاندارد ۷۰ در بین روزهای دارای انحراف مثبت برای YA‏ اکتبر ۲۰۱۷ یعنی 
۵ روز پیش از زلزله بوده است. این ناهنجاری در هر میانگین متحرکی از چند تا چند ده روزه آشکار است. 

(Akhoondzadeh, De Santis, Marchetti, Piscini, Jin, 2019)‏ نیز با بررسی و مقايسة داده‌های 555 43 
چهار ماه پیش از زلزلة ازگله با داده‌های تاریخی بلندمدت چهل‌ساله در مدل جهانی جو ECMWF‏ دریافتند 
که ۱۵ روز پیش از زلزله» ناهنجاری شدید مث مثبتی با شدت بیش از ۱۰ درجۀ سلسیوس در دمای سطح زمین 
در محدوده Y‏ درجه‌ای اطراف رومرکز این زلزله رخ داده است.(2022 (ing & Singh,‏ هم با کمک ۱۸ سال 
dale‏ سنجنده AIRS‏ از سال ۲۰۰۳ تا ۲۰۲۰ و بررسی دمای سطح زمین در منطقۀ رومرکز JJ;‏ ازگله به طول 
و عرض جغرافیایی به ترئیب Y‏ و ۱ درجه دریافتند که این پارامتر در YA‏ اکتبر و ۷ نوامبر ۲۰۱۷ دارای 


ناهنجاری مثبت بوده Al‏ 


برای درک توزیع مکانی این ناهنجاری حرارتی. اختلاف LST‏ شبگاه در روز ۲۸ اکتبر ۲۰۱۸ نسبت به 
روز قبل از خود محاسبه و به نقشه درآمد (شکل (A‏ طبق این نقشه مرز حشکی imio‏ عربستان و اوراسیا 
در غرب ایران گرم است. وجود ناحیه‌ای داغ در مرکز عراق که راستایش به رومرکز JJ;‏ ازگله می‌رسد نیز 
جالب است. گویا شدیدترین فشار صفحه‌ای از آن ناحیه وارد می‌شود. داغ‌ترین Abd‏ نقشه با ۲۶/۲ درجة 


سلسیوس اختلاف das‏ در ۲۵۰ کیلومتری جنوب غرب رومرکز در شهر فلوجة عراق واقع شده است. آنجا؛ 


داغ‌ترین ناحیه در آن روز در نیمکرۀ شمالی است و بروز ol‏ با این شدت در این منطقه از سه سال پیش از 
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شکل -٩‏ نقشة اختلاف LST‏ شبگاه در ۲۸ اکتبر ۲۰۱۷ نسبت به روز قبل از خود. 


۲۱۵۰۵۰ Map of the nighttime LST difference on October 28, 2017 compared to the 
previous day 


بررسی poled‏ ماهواره‌ای در زمان بروز این ناهنجاری, پدیدۀ ویژه‌ای را OUS‏ نمی‌دهد؛ ولی از چند 
ساعت we‏ طوفان گردوغبار سنگینی از شمال عربستان برخاسته و عراق را می‌پوشاند. افزايش ناگهانی AT‏ 
و LST‏ پیش از طوفان گردوغبار iile‏ گزارش دارد (Kayetha, Kumar, Prasad, Cervone, & Singh,‏ 
)2007 

(Wu, Zhou, Miao, & Qin, 2018)‏ نیز ناهنجاری شدیدی در شار گرمای نھان سطح (SLHF)‏ در 
غرب رومرکز و دو هفته پیش از 35 3 ازگله یافته‌اند که روز بعد با برخاستن طوفان سهمگینی از عربستان و 
عراق LLL‏ شده و دو روز بعد» محتوای (AOD) polga‏ به بیشینه رسیده است.(2022 (Jing & Singh,‏ 
هم وجود ناحیه‌ای داغ با LST CYL poli‏ در ۲۷ و ۲٩‏ اکتبر را در شمال عربستان گزارش کرده‌اند که تا 
آخر ماه میلادی به سمت شمال اپران حرکت کرده است. ایشان به رخداد طوفان در منطقه و بیشینه OAS‏ 
AOD‏ هم اشاره داشته‌اند. 


! Surface Latent Heat Flux 
? Aerosol Optical Depth 


سال سیزدهم حیدری و همکاران» بررسی پیش‌نشانگر اهنجاری حرارتی در SEE‏ ۱۳۹۹ از AS‏ کرمانشاه ۱۰۷ 


پررسی ناهنجاری حرارتی در زلزله‌های قوی ایران 


نخست. زازله‌های دستگاهی و قوی ایران از کاتالوگ USGS‏ استخراج شد (شکل ۱۰). طبق این 
فهرست. در دوره دستگاهی. YN‏ ;43 قوی L‏ بزرگای دست کم ol Jl 22 "Vo‏ دح داده Fu‏ چهار ;43 P‏ 


این فهرست در دوره عملیاتی سنجنده MODIS‏ رخ داده و می‌توان يافتەھا را روی آن‌ها نیز آزمود. 
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شکل ۱۰- نقشة زلزله‌های دورۀ دستگاهی ایران با بزرگای دست کم ۵ از کاتالوگ USGS‏ 


Fig.10. Earthquakes map of Iran in instrumental era with a magnitude of at least 6.5 from 
the USGS catalog 


برای هر یک از این چهار 35 اخیر. محدوده‌ای مربعی به طول ۲۵ کیلومتر به مرکزیت رومرکز آن‌ها 
zb‏ کاتالوگ 5 انتخاب شد و داده‌های LST‏ شبگاە از همان محصول MODIS‏ برای حدود یک سال 
پیش از هرکدام استخراج گردید و با همان dal à‏ کشف ناهنجاری حرارتی روی آن‌ها آزموده شد. نتایج 
نشان داد زلزلۀ VT‏ بم کرمان به سال ۲۰۰۳ (شکل ۰۱۱ 3335 ۷۷ محمدآباد ریگان کرمان به سال ۲۰۱۰ 


(شکل (VY‏ و البته زلزلة ازگله (شکل ۱۳) دارای پیش‌نشانگر حرارتی بوده‌اند. 


۱۰۸ نشریه Lal ae‏ و مخاطرات محیطی. جلد VP‏ شماره ۲. تابستان ۱۶۰۳ 


A: Daily Nighttime LST Averages in 2003 Bam Epicentral Region (Cloudiness <50%) 
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B: Daily Anomaly in Land Surface Temperature (LST) Series 
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شکل ۱۱- سری زمانی LST‏ شبگاه در 335 بم: -A‏ داده‌های روزانه؛ -B‏ ناهنجاری روزانه. 


Fig.11. Nighttime LST time series in Bam earthquake: A- Daily data; B- Daily anomaly 


A: Daily Nighttime LST Averages in 2010 M.A. Rigan Epicentral Region (Cloudiness <50%) 
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B: Daily Anomaly in Land Surface Temperature (LST) Series 
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شکل ۱۲- سری زمانی LST‏ شبگاه در 335 محمدآباد ریگان: -A‏ داده‌های روزانه؛ -B‏ ناهنجاری روزانه. 
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Fig.12. Nighttime LST time series in Mohammadabad Rigan earthquake: A- Daily data; 
B- Daily anomaly 


سال سیزدهم حیدری و همکاران» بررسی مر تا گر اهنجاری حرارتی در SERVE,‏ ۱۳۹۹ از AS‏ کرمانشاه ۱۰۹ 


A: Daily Nighttime LST Averages in 2017 Azgeleh Epicentral Region (Cloudiness <50%) 
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شکل ۱۳- سری 565( LST‏ شبگاه در زلزلة ازگله: -A‏ داده‌های روزانه؛ -B‏ ناهنحاری روزانه. 


Fig.13. Nighttime LST time series in Azgeleh earthquake: A- Daily data; B- Daily anomaly 


Gb‏ نمودارهای AL‏ شدیدترین ناهنجاری حرارتی در رومرکز و در سال منتهی به این سه al sly‏ حداکثر 
۵ روز پیش از وقوع آن‌ها رخ داده که بی‌سابقه‌ترین برای زلزلة ازگله است (جدول Y‏ پیش‌نشانگر حرارتی» 
تنها در زلزلة WV‏ سراوان سیستان و بلوچستان به سال ۲۰۱۳ بروز نکرده است. این زلزله. زمستانی نیست و 
عمق زیادی دارد. گفته شده که ناهنجاری حرارتی در زلزله‌های با عمق بیشتر از 1۰ کیلومتر بروز سطحی 
ندارد(2002 (Tronin, Hayakawa, & Molchanov,‏ اگرچه چندین مطالعه با روش آماری (Khalili, Alavi‏ 
L Panah, Abdollahi Eskandar, 2019)‏ به‌کارگیری شبکهة عصبی (2020 (Nekoee & Shah Hosseini,‏ بروز 


جدول ۳- نتایج آزمون کشف ناهنجاری حرارتی روی چهار UI‏ قوی و اخیر ایران. 


Table3- The results of thermal anomaly detection test on four strong and recent 
earthquakes in Iran 


Bay 1 ۳‏ عمق نمره روز : 
زلز تارب" مرجع موافق 
€ مه e we 0 km‏ 


(Choudhury, et al., 2006) à TY 13 VA Ye YANAAS بم‎ 
محمدآباد‎ 
(Akhoondzadeh, et al., 2019) A is N VAS ۳ 
ریگان‎ 
= (41) 8/0 As V۷ ANELAT سراوان‎ 


(Saber Mahani, et al., 2017) ۱۵ 0/۲ ۱۹ VY YAVA VAY از گله‎ 


۱۶۰۳ تابستان‎ CY شماره‎ OY نشریه جغرافیا و مخاطرات محیطیء جلد‎ n. 


i‏ نتیجه گیری 


امروز. افسانة باستانی گرم شدن هوا و زمین پیش از زلزله‌هاء با نام ناهنجاری حرارتی در صدها پژوهش 
به اثبات رسیده است؛ اما استفاده از روش‌های متعدد در Ges‏ اهنجاری و بررسی بازه‌های زمانی و 
محدوده‌های مکانی گوناگون. مانع از جمع‌بندی نهایی و UI‏ سامانه‌ای برای پیش‌بینی زلزله شده است. حال» 
daw 3‏ علوم داده‌کاوی در داده‌های حجیم و هوش مصنوعی. امیدها را برای راه‌اندازی چنین سامانه‌ای زنده 
کرده است. 

در این پژوهش. وجود ناهنجاری حرارتی در 3335 ازگله در سری داده‌های هواشناسی و ماهواره‌ای از سه 
سال پیش از زلزله با روش dole‏ انحراف از میانگین متحرک موردبررسی قرار گرفت و دیده شد که یکی دو 
هفته پیش از زلزله» هم در پارامترهای هواشناسی و هم در LST‏ ناهنجاری روزانه‌ای بروز می‌کند. 

بررسی آمار ایستگاه همدید تازه‌آباده OUS‏ داد که پارامترهای دمای dga‏ دمای خاک و رطوبت هواء یکی 
دو هفته مانده به زلزله ناهنجار شده‌اند؛ مثلاً دمای کمینه در VE‏ روز, دمای سطح خاک ساعت T‏ جهانی در 
۳ روز و فشار بخار ساعت ٩‏ جهانی در V‏ روز پیش از زلزله» انحراف مثبت بی‌سابقه یا کم‌سابقه‌ای داشته‌اند. 


بررسی زمانی روند LST‏ فاش کرد که بیشترین انحراف» ۱۵ روز پیش از زلزله رخ داده است. بررسی 
مکانی این ناهنجاری» روشن ساخت که در آن روز. مرز خشکی دو صفحة عربستان و اوراسیا گرم بوده و 
ناحیه‌ای Flo‏ در مرکز GLE‏ و در راستای رومرکز زلزله پدیدار شده است. از چند روز بعد» پیش‌نشانگر ابر 
زلزله )2024 (Heidari, Mazidi, & Rousta,‏ و ناهنجاری یونوسفری(2020 (Senturk, Inyurt, Sertcelik,‏ 
هم بروز نموده‌اند. 

همچنین» آزمون این فرآیند کشف ناهنجاری حرارتی روی سه زلزلۀ قوی و اخیر دیگر olal‏ نشان داد که 
دو زلزلۀ بم و محمدآباد ریگان نیز دارای ناهنجاری بی‌سابقه‌ای به ترتیب در ٦‏ و ٩‏ روز پیش از زلزله بوده‌اند. 
در پایان» تکرار چنین پژوهش‌هایی با روش‌های ماشینی کشف ناهنجاری پيشنهاد می‌شود؛ تا درنهایت بتوان با 


پایش خودکار LST‏ و L‏ پارامترهای هواشناسی همچون دمای هواء به امکان صدور هشدار زلزله رسید. 


سپاسگزاری 
از پشتیبانی دانشگاه یزد در به انجام رسیدن پایان‌نامه مطالعه پیش‌نشانگرهای زلزله ۱۳۹۹ ازگله با کمک 


سنجش ازدور, که این پژوهش, بخشی از یافته‌های OF‏ بود. صمیمانه سپاسگزاری می گردد. 


سال سیزدهم حیدری و همکاران» بررسی پیش نشانگر ناهنجاری حرارتی در زمین لرزهٌ ۱۳۹۹ از AS‏ کرمانشاه ۱۱۱ 
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